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Executive  Summary; 

The  environment  significantly  affects  minehunting  system  performance.  This  is 
especially  true  in  difficult  minehunting  conditions  to  detect  low  contrast  mines 
in  cluttered  environments.  Inexact  performance  predictions  result  in  poor  eval¬ 
uations  of  mine  clearance  which  in  turn  lead  to  inaccurate  assessments  of  re¬ 
maining  risk  to  follow-on  forces  in  the  channel  or  area.  In  this  context,  the 
present  report  provides  a  quantification  and  understanding  of  environmental 
phenomena  critical  to  acoustic  backscattering  at  minehunting  frequency.  The 
following  study  demonstrates  the  unambiguous  importance  of  seabed  interface 
roughness,  internal  bioturbation  and  upper  sediment  structure  for  object  de¬ 
tection.  In  some  cases  (e.g.  in  the  presence  of  strong  seabed  inhomogeneities), 
minehunting  sonar  detection  performances  are  believed  to  be  strongly  reduced. 
In  order  to  remedy  to  these  environmental  limitation,  it  is  clear  that  minehunt¬ 
ing  operations  need  to  quantitatively  account  for  the  environmental  conditions 
and  adapt  to  them  in  near  real  time.  The  development  of  minehunting  systems 
having  an  adaptive  behavior  as  a  function  of  the  environment  (i.e.  in  term  of 
frequency  sweep,  sensing  geometry,  multi  aspect  coverage,  etc...)  is  currently 
the  object  of  ongoing  research  at  SACLANTCEN. 
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E.  Pouliquen,  A.  P.  Lyons,  N.  G.  Pace, 
E.  Michelozzi,  L.  Muzi 


Abstract: 

Recent  backscattering  measurements  made  in  the  Gulf  of  La  Spezia  (Italy)  using 
a  sonar  operating  at  140  kHz  combined  with  thorough  seabed  interface  and  vol¬ 
ume  ground  truth  illustrate  the  dominance  of  seabed  volume  scattering*  Three- 
dimensional  fluctuation  statistics  of  density  variability  and  vertical  density  gradients, 
both  of  which  relate  directly  to  the  level  of  bioturbation  (e.g.  sea  shell  fragments, 
burrows,  pockets  of  water)  have  been  quantified  using  X-Ray  computed  tomogra¬ 
phy.  Two-dimensional  interface  roughness  spectra  have  also  been  determined  using 
a  digital  stereo  photogrammetry  system.  The  combined  ground  truth  has  allowed 
a  backscattering  model  to  be  fully  constrained.  Measured  backscattering  strength 
versus  angle  is  compared  to  a  model  that  includes  the  effects  of  varying  density 
and  sound  speed.  Data-model  comparisons  show  that  scattering  from  the  volume 
of  strongly  inhomogeneous  sediments  can  often  be  a  primary  contributor  to  seafloor 
scattering  away  from  normal  incidence. 

Resume: 

De  recentes  mesures  effectuees  dans  le  Golfe  de  La  Spezia  (Italie)  par  un  sonar  operant 
a  une  frequence  de  140  kHz  combin^es  a  une  mesure  approfondie  des  proprietes 
d’interface  et  de  volume  ont  permis  d’illustrer  la  dominance  de  la  diffusion  de  vol¬ 
ume.  Des  mesures  tri-dimensionnelles  des  fluctuations  statistiques  de  la  variation  de 
la  densite  ainsi  que  des  mesures  de  gradients  verticaux  (ces  phenomenes  etant  tous 
deux  lies  k  la  presence  de  fragments  de  coquilles,  de  terriers  et  de  poches  d’eau)  ont 
permis  de  quantifier  de  nombreux  parametres  de  volume  grace  a  une  technique  de 
tomographic  utilisant  des  rayons  X.  La  rugosite  bi-dimensionnelle  d’interface  a  ete 
determinee  grace  a  un  systeme  ster^o-photogramm^trique  numerique.  La  combinai- 
son  de  ces  deux  mesures  in-situ  a  permis  de  fournir  tous  les  parametres  necessaires 
k  un  mod^e  simple  de  retrodiffusion.  Les  mesures  angulaires  de  section  efficace  de 
r^trodiffusion  sont  compar^es  au  modele  incluant  les  effets  de  variation  de  densite 
et  de  celerite.  La  comparaison  modMe-donnees  montre  que  la  diffusion  de  volume 
de  sediments  tres  inhomogenes  pent  tres  souvent  etre  le  premier  contributeur  a  la 
diffusion  totale,  principalement  aux  incidences  non  verticales. 

Keywords:  high-frequency  acoustics,  acoustic  scattering,  seafloor  scattering. 
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1 

Introduction 


The  increasing  use  of  sonars  and  echosounders  at  very  high  frequencies  (i.e.  /  >  100 
kHz)  for  bathymetry,  seabed  characterization  and  object  detection  relies  on  the  as¬ 
sumption  that  scattering  from  the  seabed  is  essentially  caused  by  the  rough  water- 
sediment  interface  and  its  impedance  contrast.  This  assumption  is  mainly  moti¬ 
vated  by  the  (apparently)  prohibitive  high  attenuation  of  the  compressional  wave 
within  the  sediment  volume.  However,  very  high  frequency  measurements  of  seabed 
backseat tering  strength  have  not  been  able  to  validate  this  hypothesis  and  prove 
an  unambiguous  dominance  of  interface  scattering.  In  contrast  to  measurements 
below  100  kHz  that  impute  scattering  from  the  seabed  to  both  interface  and  vol¬ 
ume  scattering,  [1]-[10],  very  little  knowledge  is  available  above  100  kHz.  Stanic 
[10]  presented  backscattering  data  above  lOOkHz  from  which  it  was  impossible  to 
identify  either  a  clear  angular  or  frequency  dependence,  or  the  physical  causes  of  the 
measured  backscattering  levels. 

Until  recently,  modeling  and  validated  measurements  at  these  very  high  frequencies 
have  been  limited  by  the  difficulty  of  accurately  characterizing  seabed  properties 
at  centimeter  and  sub-centimeter  scales.  Both  interface  roughness  and  the  internal 
structure  of  the  first  tens  of  centimeters  of  the  seabed  require  accurate  measure¬ 
ment  tools  that  are  non-destructive  and  that  have  high  spatial  resolution.  These 
tools,  such  as  underwater  two-dimensional  digital  photogrammetry  [11]  and  three- 
dimensional  X-Ray  computed  tomography  [12]  [13],  have  only  become  recently  avail¬ 
able.  In  addition  to  classical  techniques  (e.g.  core  velocity  profiling  and  grain  size 
analysis),  these  new  techniques  allow  quantification  of  interface  roughness  spectra 
as  well  as  millimeter  scale  resolution  measurement  of  the  internal  sediment  density 
structure.  These  parameters  are  direct  inputs  to  acoustic  scattering  models  and 
enable  their  validation. 

Physics-based  models  designed  for  applications  below  100  kHz  have  not  yet  been 
validated  at  higher  frequencies.  The  so-called  Lambert’s  law  is  still  the  most  com¬ 
monly  used  approach  to  empirically  predict  backscattering  from  the  seabed  above 
100  kHz  [14].  This  law  implicitly  assumes  that  the  seafloor  is  a  perfectly  diffuse  re¬ 
flector  in  the  upper  space.  This  simple  approach  reduces  the  effect  of  many  possible 
seabed  scattering  mechanisms  to  a  single  and  often  physically  ambiguous  value.  For 
example,  what  appears  to  be  acoustically  bright  on  sidescan  sonar  images  could  be 
caused  by  a  seafloor  with  a  smooth  interface  and  inhomogeneous  sub-structure  in- 
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stead  of  the  expected  rough  and  hard  surface.  With  the  advent  of  new  tools  capable 
of  finely  characterizing  the  seabed  on  a  sub- centimeter  scale,  it  now  seems  natural 
to  treat  above  100  kHz  both  interface  and  volume  scattering  mechanisms  (instead  of 
simply  treating  seabed  scattering  as  a  Lambertian  process)  as  each  may  be  caused 
by  different  and  sometimes  independent  seabed  physical  processes. 

Vertical  velocity,  demsity  and  attenuation  gradients  observable  in  the  first  few  cen¬ 
timeters  of  unconsolidated  sediments  have  the  potential  to  particularly  affect  high 
frequency  scattering.  In  shallow  water  where  bioturbation  in  sediments  is  often 
present,  these  gradients  often  exist  over  vertical  distances  larger  than  the  incident 
acoustic  wavelengths,  significantly  modifying  the  local  water-sediment  impedance 
contrast  and  related  acoustic  properties.  As  a  consequence,  the  concept  of  a  re¬ 
flection  and  transmission  coefficient  between  two  homogeneous  media  both  having 
constant  average  geo-acoustic  properties  may  not  be  appropriate  at  very  high  fre¬ 
quency.  As  will  be  detailed  further,  the  presence  of  property  gradients  often  reduces 
interface  backscattering  levels  while  strongly  augmenting  the  relative  proportion  of 
volume  backscattering  to  total  backscattering. 

Although  it  seems  a  priori  counter  intuitive,  the  present  study  will  show  that,  at  very 
high  frequency,  volume  scattering  from  the  thin  upper  layer  of  unconsolidated  soft 
and  inhomogeneous  sediments  may  be  the  dominant  component  of  bottom  backscat¬ 
tering,  especially  in  the  mid-angle  range  where  dominance  of  volume  scattering  has 
often  been  observed  below  100  kHz.  The  present  work  also  demonstrates  the  im¬ 
portance  of  the  effect  of  upper  sediment  gradients  on  very  high  frequency  seabed 
backscattering.  The  next  section  describes  the  backscattering  models  used  to  ac¬ 
count  for  gradient  effects,  interface  and  volume  backscattering.  Section  3  presents 
the  fine  scale  gcoacoustic  characterization  of  four  sites  in  the  Gulf  of  La  Spezia  (Italy) 
where  density  gradients  and  density  fluctuations  were  observed  and  quantified.  In 
section  4,  a  comparison  between  the  measured  and  predicted  backscattering  strength 
(BS)  at  these  four  sites  is  made.  Section  5  presents  a  summary  and  a  discussion  on 
the  obtained  results. 
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Model  description 


To  predict  the  relative  proportion  of  sediment  volume  backscattering  against  in¬ 
terface  backscattering  above  100  kHz,  several  approaches  are  available.  Classical 
approaches  like  the  ones  presented  by  Jackson  [15]  and  Yamamoto  [16]  have  been 
retained  in  this  work  mainly  for  the  sake  of  simplicity.  These  relatively  simple  models 
have  been  validated  below  100  kHz  and  are  based  on  the  tangent  plane  approxima¬ 
tion  (TPA)  and  on  the  small  perturbation  theory  (PT)  to  treat  backscattering  from 
the  water-sediment  interface  [15]  and  on  the  small  perturbation  model  to  predict  vol¬ 
ume  backscattering  [16]-  [19].  The  TPA  is  an  intuitive  high  frequency  approximation 
derived  from  physical  optics.  It  becomes  inaccurate  when  shadowing  occurs  (i.e  typ¬ 
ically  when  the  grazing  angle  becomes  small).  Regarding  the  perturbation  theory  as 
applied  to  interface  scattering,  it  can  be  seen  as  a  lower  frequency,  low  grazing  angle 
approach.  The  conditions  for  which  each  of  these  theories  break  down  for  pressure 
release  surfaces  are  discussed  in  detail  in  [20]  [21]  for  Gaussian  roughness  spectra. 


2.1  Interface  scattering 


The  backscattering  strength  is  assumed  to  be  the  incoherent  sum  of  the  interface  and 
volume  backscattering  strengths.  Each  interface  and  volume  component  is  therefore 
treated  separately.  It  is  also  assumed  that  the  coherent  component  of  the  backscat- 
tered  field  is  zero. 


2.1.1  Tangent  plane  approximation 

To  model  the  interface  backscattering  cross  section  for  the  case  of  no  gradient,  we  use 
the  classical  tangent  plane  approximation  (or  Kirchhoff  approximation).  Following 
Jackson  [15]  but  using  the  reflection  coeflicient  evaluated  at  the  grazing  angle  rather 
than  at  normal  incidence,  the  TPA  yields  for  the  interface  scattering  cross  section: 
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where  ka  is  the  wavenumber  in  the  water,  6  is  the  grazing  angle,  and  R\2{0)  is 
the  plane  wave  reflection  coefficient  between  the  water  of  impedance  and  the 
sediment  of  impedance  psCg  (in  the  case  of  no  gradient).  In  this  study,  for  the  sake  of 
simplicity,  we  assume  an  isotropic  power-law  spectrum  for  roughness  of  the  following 
form: 


W{k)  =  (2) 

As  described  in  greater  detail  in  [15],  the  structure  function  D  present  in  Eq.  1  and 
the  roughness  power  spectrum  arc  connected  by  the  following  transformation: 


/  +  00  /•+00 

/  (1  —  cos(k.r))W'(k)d^fc. 

-00  J—oo 


(3) 


In  the  case  of  the  power  law  spectnim  of  Eq.2,  the  structure  function  takes  the 
simple  form: 


D{v)=dir-<-\ 

with 


(4) 


Cl  =  I27r/5r(3  -  7/2)2-(^-2)]/[(7/2  -  1)(2  -  7/2)r(7/2)].  (5) 

The  technique  used  to  characterize  the  roughness  power  spectrum  and  obtain  values 
for  P  and  7  is  described  in  section  3.3. 


2.1.2  Perturbation  theory 

When  the  small  scale  RMS  relief  is  much  smaller  than  the  acoustic  wavelength,  the 
Rayleigh-Rice  perturbation  approximation  (or  simply  the  perturbation  approxima¬ 
tion)  can  be  applied  to  the  small  scale  portion  of  the  water-sediment  interface  [22]. 
In  this  paper,  the  backscattering  strength  using  the  first  order  perturbation  theory 
also  follows  Jackson’s  formulation  for  a  penetrable  two-fluid  interface  [15]  and  is 
expressed  as 


(y,{9)  =  sm\0)Fie,  i/,m)W"(2fc„cos(0)), 


(6) 
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with, 


F(0,  z/,  m) 


'  f(m  — 1)^  COS^(0) +771^—1^ 

[msin(^)  +  (t^-2-cos2(^))V2]4 

< 

[(m  — 1)^  cos^ 

k.  [(1— m2)  cos2(0)+m2— l/~2^j2 


e>ec 

9  <00. 


(7) 


The  parameter  m  in  this  equation  is  the  sediment-water  density  ratio  (m  ps/pw)^^ 
is  the  velocity  ratio  {u  =  Cs/c^)  and  9c  the  critical  angle.  The  expression  in  Eq.  6 
is  also  referred  to  as  the  small  scale  backscattering  cross  section.  Assuming  that  the 
slope  of  the  large-scale  surface  is  a  Gaussian-distributed  random  variable,  Jackson 
averages  this  cross  section  to  account  for  shadowing  [15].  This  assumes  a  cutoff 
wavenumber  in  the  roughness  spectrum  that  is  somewhat  arbitrary  but  that  has 
the  merit  of  providing  a  solution  that  accounts  for  the  shadowing  effect.  Following 
Eqs.  19-24  from  [15],  one  obtains  a  low  grazing  angle  solution  to  the  backscattering 
cross  section,  setting  the  RMS  slope  of  the  large-scale  surface  arbitrarily  to  0.08 
radian. 


2.2  Volume  scattering 

To  treat  scattering  from  sediment  density  and  velocity  fluctuations,  the  expression 
for  the  volume  backscattering  strength  (i.e.  the  differential  backscattering  cross 
section  per  unit  surface)  is  taken  from  Yamamoto  [16]:  . 


cos(0') 

4q' 


sin^(0')  cos^(0) 
sin^(0)  cos^(^') 


(8) 


The  quantity  a  stands  for  the  sediment  compressional  attenuation,  T12  is  the  water- 
sediment  transmission  coeflftcient,  r2i  is  the  sediment- water  transmission  coefficient, 
9  is  the  incident  angle,  0'  is  the  refracted  grazing  angle.  This  expression  is  obtained 
for  a  smooth  water-sediment  interface.  Ideally,  the  scattering  effect  at  the  crossing 
of  the  rough  interface  should  be  taken  into  account  as  is  done  by  Ivakin  [19]  but 
for  smooth  surfaces,  it  has  a  negligible  effect  on  the  resulting  volume  backscattering 
strength.  The  scattered  differential  cross  section  per  unit  volume  cr^  is  obtained  from 
Chernov  [17]  and  is  based  on  the  first  order  perturbation  theory.  With  the  assump¬ 
tion  that  relative  density  variations  and  relative  velocity  variations  are  correlated 
accordingly. 


Sc  Sp 

—  =  r/— , 

Cs  ps 


(9) 


5- 


SaCLANTCEN  SR-342 


the  local  volume  baekscattering  cross  section  <7„  can  be  expressed  as 


a^{e)  =  2nkt{l  +  n?Sp{2ka),  (10) 

where  Sp{2ka)  is  the  three-dimensional  (3-D)  power  spectrum  of  the  relative  density 
fluctuations  evaluated  at  the  Bragg  wavenumber,  as  was  done  for  interface  scatter¬ 
ing.  The  value  of  Sp{2ka)  is  estimated  from  X-Ray  computed  tomography  scans  as 
described  below.  Also,  as  only  density  fluctuation  measurements  are  available,  the 
value  of  r]  needs  to  be  flxed  arbitrarily.  It  is  fixed  to  1/3  in  accordance  with  average 
values  measured  by  Yamamoto  [16]  [23]  for  unconsolidated  sediments. 


2.3  Effect  of  vertical  gradients 

As  the  acoustic  frequency  increases  (and  wavelength  decreases),  it  is  no  longer  pos¬ 
sible  to  consider  the  area  between  water  and  sediment  as  a  well-marked  interface 
between  two  media.  At  very  high  frequency,  vertical  gradients  may  be  of  larger 
scales  than  the  actual  acoustic  wavelength  and  may  significantly  affect  reflection 
and  transmission  from  and  through  such  a  medium.  Here,  the  effect  of  gradient 
is  considered  as  only  affecting  the  reflection  and  transmission  coefficient.  In  the 
presence  of  gradient,  an  effective  reflection  coefficient  is  substituted  to  R12 

in  Eq.l.  The  effective  transmission  coefficients  and  replace  T12  and  T21 
respectively  in  Eq.8.  In  Eq.7,  the  density  ratio  m  and  velocity  ratios  are  replaced 
by  effective  ratios,  rricff  and  Uejf  respectively. 

There  are  several  approaches  for  computing  Rl{^  and  This  can  be  done 

following  Brekhovskikh’s  approximation  of  the  reflection  coefficient  for  a  gradient 
described  by  an  arbitrary  law  [6]  and  the  particular  case  of  a  density  profile  of  the 
following  form  (see  Lyons  [25]): 


p{^)  =  Pw  if  2:  <  0  (11) 

P(^)=  if^>0  (12) 

However,  in  the  presence  of  both  density  and  velocity  gradients,  a  numerical  com¬ 
putation  is  at  present  the  only  tractable  way  to  obtain  these  coefficients.  In  the 
present  case,  velocity  is  assumed  to  behave  like  density  (p  in  Eq.l2  is  replaced  by 
c).  The  OASES  reflection  coefficient  module  (OASR  [26])  is  used  to  compute  the 
plane  wave  coefficients.  It  should  be  noted  that,  in  the  case  of  gradient,  the  scatter¬ 
ing  cross  section  obtained  with  the  TPA  and  that  obtained  with  small  perturbation 
theory  are  expected  to  differ  because  of  a  difference  in  computing  R\^/  and  the  term 
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F(05  i^e//5^e//)  ill  EQ-  7.  In  the  former  case,  Ri2^  is  obtained  from  OASR  and  in 
the  latter  z/e//5^e//)  is  obtained  using  rrieff  =  po/pw  and  Ue/j  =  cq/c^e;.  The 
following  section  illustrates  the  effect  of  gradient  on  the  reflection  coefficient  for  the 
considered  study  sites. 
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3 

Seabed  geo-acoustic  properties 


Four  sites  in  the  Gulf  of  La  Spozia  (Italy)  have  been  considered  in  this  study:  Tellaro 
(TE),  Vcnere  Azzura  (VA),  Porto  Vcnere  (PV)  and  Punta  della  Mariella  (PdM).  The 
sites  have  a  water  depth  which  varies  from  about  10  to  20  m.  For  each  of  the  sites, 
diver  cores  and  stereo-photographs  were  obtained.  At  two  of  the  locations,  PV  and 
PdM,  visibility  was  too  poor  to  obtain  good  quality  stereo-images  of  the  seabed 
interface.  The  stereo  photographs  were  further  processed  to  obtain  two-dimensional 
height  fields  and  the  cores  were  X-ray  CT  scanned  (Fig.  1)). 


3.1  Core  data  description 

The  four  samples  that  were  collected  have  lower  average  grain  size  than  those  pre¬ 
sented  in [13],  ranging  from  very  fine  sand  to  silty-clay.  The  grain  size  analysis 
presented  in  Table  1  reveals  for  each  of  them  very  different  combinations  of  sand, 
silt  and  clay.  To  reduce  disturbance,  the  cores  were  taken  manually  by  divers.  All 
were  taken  with  the  longitudinal  axis  of  the  core  liner  kept  horizontal  (hereafter 
referenced  as  the  Z-axis).  This  uncommon  geometry  was  justified  by  the  difficulty 
of  scanning  a  vertical  core.  Since  the  internal  diameter  of  the  gantry  around  which 
X-Ray  emmiters  and  detectors  rotate  is  fixed,  a  vertical  core  with  the  requirement 
that  its  original  in  situ  orientation  be  kept  would  have  had  to  be  short  enough  to 
pass  through.  As  this  study  concerns  very  high  frequency  volume  scattering  (i.e. 
typically  above  100  kHz,  corresponding  to  centimetric  to  sub-centimetric  acoustic 
wavelengths),  a  core  diameter  of  10  cm  was  found  to  be  sufficient  to  capture  sed¬ 
iments  over  several  wavelengths  in  the  vortical  direction  (i.e.  the  Y-axis).  The 
horizontal  sampling  also  enabled  an  analysis  of  the  horizontal  stationarity  of  density 
properties. 

X-ray  CT  scans  were  performed  within  hours  after  in-situ  sampling  to  reduce  com¬ 
paction  and  decomposition.  Although  the  cores  were  kept  unopened  during  the 
entire  procedure  there  exists  the  possibility  that  some  of  the  cores  may  not  have  re¬ 
mained  at  ambient  pressure  until  the  CT  scans  took  place.  The  cores  were  scanned 
using  a  Picker  PQ2000S  scanner.  Voxel  dimension  varied  between  cores  as  devseribed 
in  Table  2. 

Although  different  types  of  interaction  are  possible,  the  relatively  low  energy  of 
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Figure  1  Illustration  of  the  procedure  employed  to  perform  X-Ray  tomography  on 
the  horizontal  cores. 


Table  1  Grain  size  properties  at  the  four  study  sites 


Water  depth 
(m) 

Sand/Silt /Clay 

%/%/% 

Mean  grain  diameter 
(mm) 

TE 

14.2 

75/15/10 

0.0507 

VA 

9.8 

20/70/10 

0.0145 

PV 

12.3 

00/40/60 

0.0018 

PdM 

21.5 

15/50/35 

0.0068 

X-rays  used  by  CT  scanning  reduces  the  number  of  relevant  phenomena  to  photo¬ 
electric  effect  and  Compton  scattering.  It  has  been  determined  that  for  seafloor 
sediments,  there  is  a  linear  dependence  between  bulk  density  and  CT  scan  numbers 
(in  Hounsfield  units)  [12].  The  calibration  approach  discussed  in  [12]  is  used  to 
determine  the  relationship  between  Hounsfield  units  and  bulk  density. 

Subsequent  data  processing  similar  to  that  described  by  Tang  [13]  was  carried  out 
to  obtain  estimates  of  power  spectra.  As  described  in  detail  in  [13],  the  three- 
dimensional  density  power  spectrum  is  effectively  computed  as  follows: 


AxAyAz 

{27r)m^NyN, 

Na:-l,Ny-l,N,-l 


M\ 


E 

/,m,n=0 


Hi 


p{l,m,n)  -  p 


(13) 
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Table  2  Geometrical  properties  of  each  X-ray  CT  data  set. 


Voxel  size  (mm^) 
(mm^) 

Ax  X  Ay  X  Az 

Slice  pixel 
number 

Nx  X  Ny 

Number 

of  slices 
N, 

TE 

.35  X  .35  X  2 

290  X  290 

293 

VA 

.25  X  .25  X  2 

318  X  318 

137 

PV 

•35  X  .35  X  2 

318  X  318 

119 

PdM 

.25  X  .25  X  2 

400  X  400 

183 

^-i2Tr[lu/Nx+rnv/Ny^nw/Nz]  j2^ 

The  term  p  is  the  average  density  within  the  core  sample,  away  from  the  edges 
where  compaction  may  occur.  A  three  dimensional  window  function  Ht^rn,n  is  used 
to  reduce  leakage  in  the  spectrum  estimation.  As  in  Tang  [13],  we  use  a  three 
dimensional  Hanning  function  H{P). 

H{P)  =  !-((//-  1.)  cos(jrP)  +  p)  (14) 

where  P  is  a  scaled  radius  which  does  not  exceed  unity. 

The  term  fi  is  fixed  to  0.5,  H  takes  the  form: 

H{P)  =  i(l  +  cos(7rP))  (16) 

The  window  function  is  normalized  such  that 

Tinr  I  /  1-  (17) 

The  description  of  the  internal  structure  of  sediments  using  data  derived  from  X-ray 
CT  techniques  has  been  the  object  of  recent  studies  (mainly  by  Orsi[12]  and  more 
recently  by  [13]).  In  [13],  the  analysis  of  the  density  spectra  of  sandy  sediments 
showed  three  types  of  variability:  the  intrinsic  variability  of  the  sand  matrix,  shells 
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and  mud  mixtures.  The  three  dimensional  power  spectra  of  density  fluctuation  were 
well  modeled  by  a  two  component  power  law,  and  in  contrast  to  previous  assump¬ 
tions,  an  exponential  correlation  function  was  found  not  to  be  suitable  for  these 
sediments.  Even  though  these  conclusions  might  be  valid  for  some  sandy  sediments, 
a  number  of  questions  still  remain  open.  Despite  the  apparent  isotropic  nature  of 
the  sediments  studied  by  Tang  [13],  organized  distributions  of  individual  scatters 
(shells  or  mud  intrusion)  or  the  presence  of  layering  can  often  be  observed.  Their 
effects  on  acoustic  scattering  are  not  yet  understood.  While  modified  power  law 
spectra  are  possibly  suitable  for  sandy  sediments,  it  is  still  of  interest  to  determine 
whether  or  not  this  spectral  description  is  also  suitable  for  other  sediment  types 
especially  to  those  with  finer  grain  sizes  which  are  more  likely  to  be  affected  by  a 
strong  heterogeneity  and/or  bioturbation. 

There  is  a  limitation  in  using  a  single  core  sample  to  characterize  the  spectral  nature 
of  highly  bioturbated  sediments,  as  is  the  case  in  this  study.  Even  though  low  spatial 
frequencies  are  represented  in  one  core,  information  from  a  single  data  set  may  not 
be  sufficient  to  obtain  a  well-determined  average  power  spectrum.  In  other  words, 
the  lower  the  spatial  frequency,  the  smaller  the  relative  statistical  population  and 
therefore,  the  less  well  defined  the  average  power  spectrum  for  larger-scale  fluctua¬ 
tions.  To  ensure  proper  ensemble  averages,  larger  data  sets  in  size  and  number  from 
a  statistically  homogeneous  region  are  required  but  are  not  available  in  the  present 
study.  Also,  the  Gaussian  nature  of  the  density  distribution  may  be  questioned  as 
mud  intrusion  and  shell  pieces  may  not  be  considered  as  a  simple  additive  process 
to  the  intrinsic  mixture  of  mineral  grains  and  water. 

In  the  presence  of  a  three-dimensional  isotropy  of  the  density  field,  the  spectrum 
Sp  can  be  radially  averaged  to  obtain  a  smoother  spectrum  Sp.  By  switching  to 
spherical  coordinates,  the  spectrum  Sp  averaged  over  a  wavenumber  width  of  Ak  is 
obtained  as  follows: 


r27r  r7r/2  rk-\-Ak 

X  /  /  Sp{k^(f)^'ilj)k^dkcos{4>)d^dil; 

J <f>=—7T/2  J k 


(18) 


The  one-dimensional  density  power  spectra  5^,  Sy  respectively  in  the  X,  Y  and 
Z  direction  are  computed  by  double  integration  of  the  three  dimensional  density 
power  spectrum  Su^v,w’ 


Ny-1,N,-1 

Su  —  AkyAkz  ^  ^  0^) 

m,Ti=0 
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Sv  ~  Akx^kz 


Sjt)  =  Akx^ky 


^  ^  SujV^W 

/,n=0 

Nr.-\,Ny-l 

(20) 

l^m—0 

(21) 

For  each  core,  three  one-dimensional  spectra  and  a  radially  averaged  spectrum  are 
computed.  The  results  are  given  in  the  following  subsections  and  compiled  in  Table  3. 


3.1.1  Tellaro 

Sediments  at  this  site  arc  silty-sands  with  only  few  small  shell  fragments.  The  slice  of 
the  two  dimensional  spectrum  at  =  0  (Fig.  2)  shows  an  anisotropy.  This  deviation 
can  be  due  to  an  error  of  the  estimation  of  the  vertical  axis  or  compaction  during  the 
coring  process.  Anisotropy  is  not  observed  in  the  two  other  planes.  Away  from  the 
lower  spatial  frequency  region  of  the  spectrum  which  displays  some  fluctuation  due 
to  under-representation  of  the  statistical  population,  the  one-dimensional  marginal 
spectra  show  only  small  discrepancies  with  each  other  at  higher  spatial  frequencies 
(Fig.  3).  The  slope  of  the  spectrum  is  steeper  above  log\o{k)  =  3.5  than  the  slopes 
observed  by  Tang  for  sandy  sediments  [13],  Figure  4  shows  the  behavior  of  the 
radially  averaged  three-dimensional  spectrum.  It  clearly  does  not  behave  like  a 
modified  power  law-term  of  the  form: 


Sp\mode.l-iD{k)  -  (^2) 

Emphasis  in  the  fitting  process  was  made  on  the  higher  frequency  part  of  the 
spectrum  because  the  lower  frequency  part  of  it  suffers  from  statistical  under¬ 
representation.  A  two  component  modified  power  law,  as  suggested  by  Tang  [13], 
does  not  give  a  significantly  improved  fit  at  the  expense  of  adding  three  more 
parameters.  The  solid  line  in  Fig.  4  is  the  best  fit  to  the  measured  spectrum 
(aj  =  6.0  X  10”^m,  li  =  2mm  and  71  =  2.2)  in  a  least  square  sense.  For  the 
same  reasons,  many  other  combinations  of  analytical  forms  do  not  give  a  better  fit 
over  the  entire  spectral  band. 


3.1.2  Venera  Azzura 


Sediments  at  Venere  Azzura  are  sandy-silts,  with  the  presence  of  shell  pieces  and 
some  sand.  There  is  no  visible  presence  of  anisotropy  in  any  of  the  three  planes  of  ref¬ 
erence  (Fig.  5).  This  is  confirmed  by  the  behavior  of  the  individual  one-dimensional 
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Figure  2  Slice  of  the  three  dimensional  density  power  spectrum  Sp{kx^  ky^  kz  =  0) 
from  the  Tellaro  dataset.  Units  are  in  dB  re  1  m^. 


spectra  as  displayed  in  Fig.  6.  As  for  the  Tellaro  core  analysis,  a  drop-off  is  also 
observable  around  logio{k)  =  3.5.  In  the  lower  spatial  frequency  region,  the  slope  of 
the  power  spectrum  tends  to  zero. 

In  Fig.  7,  the  best  fit  obtained  using  the  modified  power  law  of  Eq.  22  gives  more 
satisfactory  results  than  for  Tellaro  (ai  =  6.0x  li  =  2mm  and  71  =  2.2).  The 

modified  power  law  seems  a  good  candidate  to  represent  this  sediment  structure. 


3.1.3  Porto  Venere 


Sediments  at  Porto  Venere  are  clayey-silts  with  a  strong  presence  of  shell  pieces, 
water  pockets,  worm  burrows,  and  small  stones.  Figure  8  shows  a  false  color  three- 
dimensional  reconstruction  of  the  density  field  of  the  PV  horizontal  core.  Although, 
the  upper  part  of  the  core  displays  the  water-sediment  interface,  the  spectral  analysis 
was  performed  excluding  this  region.  This  reconstruction  illustrates  the  high  spatial 
variability  of  density,  both  horizontally  and  vertically.  Figure  9  is  an  artificial  view 
of  the  internal  structure  of  the  density  field  made  by  removing  the  mid-density  values 
and  keeping  visible  the  low  and  high  density  values  (respectively  the  water  pockets 
and  the  burrows  in  red  and  the  shell  pieces  in  yellow).  This  representation  shows 
that  1)  shell  pieces  are  uniformly  distributed  in  the  core,  2)  the  low-density  regions 
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Tellaro 


Figure  3  Three  one-dimensional  spectra  5„,  from  the  Tellaro  dataset.  Units 

are  in  dB  re  1  m. 


are  not  uniformly  distributed  over  the  sediment  sample,  3)  the  shape  and  length  of 
the  low  density  regions  (most  likely  burrows)  are  not  completely  random  and  tend 
to  be  oriented  more  vertically,  and  4)  the  size  of  the  low  density  regions  is  an  order 
of  magnitude  larger  than  the  size  of  the  shell  pieces. 

The  density  power  spectrum  from  this  highly  inhomogeneous  medium  is  isotropic  as 
can  be  seen  in  Fig.  10.  In  Fig.  11,  the  three  one-dimensional  power  spectra  overlay 
each  other,  except  for  the  lowest  spatial  frequency  region  where  the  values  become 
more  random.  The  similarity  of  the  one-dimensional  spectra  for  the  three  orthogonal 
directions  is  strong  evidence  of  isotropy,  despite  the  presence  of  burrows  that  have 
an  almost  vertical  orientation.  A  drop-off  in  the  spectra  is  again  observed  around 
^ogiQ{k)  =  3.5  but  is  less  strong  than  for  the  two  previous  datasets.  The  radially 
averaged  three-dimensional  spectrum  shown  in  Fig.  12  is  well  fit  by  a  modified  power 
law  above  log\Q{k)  —  2.3.  The  model  values  are  ai  =  6.0  x  lO'^^m,  li  =  1.4mm 
and  7i  =  2.8).  The  dominance  of  the  intrinsic  density  variations  in  the  silty-clay 
sediment  can  explain  this  behavior  in  the  high  wavenumber  portion  of  the  spectrum. 
Below  log\o{k)  =  2.3,  the  large-scale  fluctuations  caused  by  the  burrows  and  shell 
pieces  seem  to  be  dominating  the  internal  density  structure.  Curiously,  oscillations 
of  the  PV  spectrum  in  the  low  spatial  frequency  part  do  not  seem  to  be  larger  in 
amplitude  than  those  observed  in  the  spectrum  of  TE  where  there  is  practically  no 
large  scale  fluctuation  (only  a  few  small  shell  pieces  and  no  significant  burrows). 
One  may  therefore  strongly  question  the  cause  of  these  oscillations  of  the  spectrum 
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Figure  4  Radially  averaged  three  dimensional  spectrum  Sp  from  the 
(dots).  The  solid  line  corresponds  to  the  best  fit  of  the  model.  Units 


Tellaro  dataset 
are  in  dB  re  1 


in  the  small  wavenumber  region.  These  oscillations  may  simply  be  due  to  the  limits 
imposed  by  the  actual  physical  size  of  the  dataset. 


3.1.4  Punta  della  Mariella 


The  sediments  at  Punta  della  Mariella  are  a  mixture  of  clay  and  silt,  with  also  a 
strong  presence  of  shell  pieces,  water  pockets  and  bioturbation.  This  is  the  most 
inhomogeneous  site  of  the  study.  The  density  fluctuation  spectrum  of  the  core  taken 
at  this  site  shows  near  perfect  isotropy.  The  one-dimensional  spectra,  taken  along 
the  three  main  axes,  all  behave  identically  as  displayed  in  Fig.  13,  A  modified  power 
law  does  not  fit  well  the  three-dimensional  radially  averaged  spectrum  in  Fig.  14, 
The  best  model  values  are  ai  =  1.0  X  10“^m,  li  =  1.2mm  and  71  =  3.0.  A  mismatch 
between  model  and  data  can  again  be  observed  in  the  low  wavenumber  portion  of 
the  spectrum,  as  seen  in  the  spectra  of  TE  and  PV. 
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15000 


Figure  5  Slice  of  the  three  dimensional  density  power  speetrum  Sp{kx,  ky,  k^  =  0) 
from  the  Venere  Azzura  dataset.  Units  are  in  dB  re  1  m^. 


3.2  Density  gradient  estimation 


Computation  of  gradients  is  a  difficult  task  because  a)  the  X-Z  plane  of  the  core  ma¬ 
trix  is  not  perfectly  horizontal  and  b)  because  of  local  roughness,  all  the  individual 
profiles  p(l,0  <  m  <  Ny,n)  for  given  combinations  of  I  and  n  may  not  be  system¬ 
atically  aligned  with  each  other.  Therefore,  the  first  step  in  gradient  computation 
constists  in  aligning  all  the  vertical  profiles  p{l,0  <  m  <  Ny,n)  with  each  other. 
For  each  individual  profile,  the  depth  where  the  density  crosses  l.lg/cm^  is  defined 
as  being  at  a  nominal  ’altitude’  of  — 0.1cm.  Once  aligned,  the  second  step  consists 
in  averaging  the  profiles  as  a  function  of  depth.  Figures  15,  16,  17  and  18  present 
the  respective  averaged  profiles  at  the  four  considered  sites.  Each  of  them  displays 
a  distinct  gradient  pattern.  The  least  pronounced  gradient  is  at  Tellaro,  and  the 
most  pronounced  one  is  at  Punta  della  Mariella.  The  fit  to  the  theoretical  profile 
proposed  in  eq.  12  is  done  manually.  Table  3  at  the  end  of  this  section  displays  the 
chosen  values  of  pQ,  ps  and  a  for  each  site. 


3.3  Interface  roughness  characterization 


Measurement  of  roughness  power  spectra  W  were  obtained  using  a  digital  stereo 
photogrammetry  system  developed  at  SACLANTCEN.  Providing  visibility  condi- 
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Figure  6  One- dimensional  spectra  Su,  Sv,  Syj  from  the  Venere  Azzura  dataset. 
Units  are  in  dB  re  1  m. 


tions  are  good,  the  system  gives  a  snapshot  estimate  of  the  two-dimensional  water- 
sediment  interface  height  field  for  a  limited  area  (typically  less  than  one  square 
meter).  The  procedure  employed  to  retrieve  roughness  spectra  is  described  in  de¬ 
tail  in  [11].  Prom  the  height  field  estimate,  the  roughness  power  spectrum  can  be 
computed  (and  subsequently  values  of  ^  and  7  can  be  estimated).  Given  the  lim¬ 
ited  number  of  measurements  per  site,  stationarity  over  a  larger  area  is  generally 
assumed.  Visibility  was  sufficient  only  at  TE  and  VA  to  retrieve  roughness  spectra 
with  a  good  degree  of  confidence.  At  the  two  other  sites  PV  and  PdM,  visibility  was 
too  poor  to  enable  a  satisfactory  stereo  reconstruction. 


Figures  19  and  20  give  a  qualitative  idea  of  interface  roughness  present  at  TE  and 
VA  respectively.  One  can  notice  the  presence  of  burrow  holes  at  VA  and  their 
absence  at  TE.  There  is  a  slight  roughness  anisotropy  due  to  a  light  ripple  field  at 
TE  but  acoustic  backscattering  at  this  site  did  not  show  any  significant  difference 
between  along  and  across  the  ripple  crests.  Figure  21  displays  the  retrieved  height 
field  at  Venere  Azzura.  The  magnitude  of  roughness  at  TE  and  VA  is  low,  less  than 
one  centimeter  RMS  over  the  reconstructed  height  field.  Assuming  isotropy  at  both 
sites,  P  and  7,  are  estimated  as  defined  by  Eq.  2.  It  is  made  by  performing  a  least 
square  fit  to  the  radially  averaged  roughness  spectrum: 
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Venere  Azzuro 


Figure  7  Radially  averaged  three  dimensional  spectrum  Sp  from  the  Venere  Azzura 
dataset  (dots).  The  solid  line  corresponds  to  the  best  fit  of  the  model.  Units  are  in 
dB  re  1  . 


W{k)  = 
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2it  rk+Ak 


kdk  Jij) 


r  f 

Jtb=0  Jk 


Sp{k,'ip)kdkd'il). 


(23) 


Figures  22  and  23  show  the  roughness  power  spectra  and  their  corresponding  best 
fit  using  a  power  law  spectrum.  The  fit  is  good  above  0.06ct/c/es/cm  (i.e.  for 
wavelengths  smaller  than  approximately  20  cm). 


3,4  Seabed  geo-acoustic  properties  summary 


By  combining  X-Ray  CT  scans,  stereo  photogrammetry  techniques  and  grain  size 
analysis,  most  geo-acoustic  sediment  properties  needed  for  input  to  the  model  pre¬ 
sented  in  section  2  have  been  characterized.  They  are  compiled  in  Table  3.  Param¬ 
eters  in  parenthesis  were  those  estimated  without  a  dedicated  measurement.  Atten¬ 
uation  in  particular  was  not  measured  but  estimated  from  grain  size  measurements 
using  Hamilton  [27].  Because  roughness  estimates  at  PV  and  PdM  were  not  avail¬ 
able  due  to  poor  visibility  conditions,  /3  and  y  were  assumed  identical  to  the  values 
measured  at  VA.  The  value  of  Sp{2ka)  is  taken  from  the  radially  averaged  spectra  Sp 
at  the  spatial  frequency  of  interest  {2ka  =  lllSrad/m  at  UOkHz).  This  ensemble 
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Figure  8  False  color  3-D  reconstruction  of  horizontal  core  at  PV;  The  size  of  the 
core  is  approximately  10x10x22  cm^. 


of  geo-acoustic  properties  enables  the  backscattering  strength  model  presented  in 
Section2  to  be  almost  fully  constrained,  requiring  only  a  few  of  the  parameters  to 
be  estimated. 
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Figure  9  3-D  visualization  of  burrows  and  water  pockets  (in  red)  and  shell  pieces 
(in  yellow)  in  the  horizontal  core  of  Porto  Venere 
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Figure  10  Slice  of  the  three  dimensional  density  power  spectrum  Sp{kx^  ky^  kz  =  0) 
from  the  Porto  Venere  dataset.  Units  are  in  dB  re  1  m^. 
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Figure  11  One  dimensional  spectra  Suf  Sy,  of  Porto  Venere  along  the  three 
axis.  Units  are  in  dB  re  1  m. 


Porto  Venere 


Figure  12  Radially  averaged  three  dimensional  spectrum  Sp  from  the  Porto  Venere 
dataset  (dots).  The  solid  line  corresponds  to  the  best  fit  of  the  model.  Units  are  in 
dB  re  1  m^. 
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Punta  Della  Manella 


Figure  13  One  dimensional  spectra  W(k)  of  Punta  della  Mariella  along  the  three 
axis.  Units  are  in  dB  re  1  m 


Punta  Della  Mariella 


Figure  14  Radially  averaged  three  dimensional  spectrum  Sp  of  Punta  della 
Mariella  (dots).  The  solid  line  corresponds  to  the  best  fit  of  the  model  spectrum 
at  high  wavenumbers.  Units  are  in  dB  re  1  m® . 
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Tellaro 


Figure  15  Average  density  profile ^  intervals  of  confidence  and  theoretical  profile 
as  a  function  of  depth  within  the  sediment  of  Tellaro 


Venere  Azzura 
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Figure  16  Average  density  profile,  intervals  of  confidence  and  theoretical  profile 
as  a  function  of  depth  within  the  sediment  of  Venere  Azzura 
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Porto  Venere 


Figure  17  Average  density  profile,  intervals  of  confidence  and  theoretical  profile 
as  a  function  of  depth  within  the  sediment  of  Porto  Venere 


Punto  Della  Mariella 


Figure  18  Average  density  profile,  intervals  of  confidence  and  theoretical  profile 
as  a  function  of  depth  within  the  sediment  of  Punta  della  Mariella 
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Figure  21  30  cm  x  SO  cm  digital  elevation  model  at  Venere  Azzura  obtained  using 
the  stereo  photo^grammetry  technique. 


Tellaro 


Figure  22  Roughness  power  spectrum  at  Tellaro,  The  solid  line  corresponds  to  the 
best  fit  with  the  radially  averaged  power  spectrum  W.  For  comparison^  the  dashed 
line  corresponds  to  the  value  of  the  2-0  roughness  spectrum  W  at  ky  =  0. 
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Venere_Az2ura 


Figure  23  Roughness  poPouldul  wer  spectrum  at  Venere  Azzura.  The  solid  line 
corresponds  to  the  best  fit  with  the  radially  averaged  power  spectrum  W ,  For  com¬ 
parison^  the  dashed  line  corresponds  to  the  value  of  the  2-D  roughness  spectrum  W 
at  ky  =  0. 


-27- 


SaCLANTCEN  SR-342 


Table  3  Geoacoustic  properties  of  the  four  study  sites.  Parameters  in  parenthesis 
are  estimates. 


TE 

VA 

PV 

PdM 

Voliune  Parameters 

Mean  grain  diameter 

(mm) 

0.0507 

0.0145 

0.0018 

0.0068 

Sand/Silt/Clay 

%/%/% 

75/15/10 

20/70/10 

00/40/60 

15/50/35 

Density  Lab.  Measurement 

(g/crri^) 

1.85 

1.65 

1.5 

1.7 

Density  X-Ray  Measurement  ps 

(g/cm^) 

1.97 

1.99 

1.61 

1.81 

Compress.  Attenuation  a 

(dB/m/kHz) 

(0.6) 

(0.04) 

(.10) 

Sound  speed  Cs 

(m/s) 

1650 

1490 

1530 

Sp{2ka)  Q140kHz 

-147.01 

bllbji 

-142.1 

-138.5 

Voliime  Scat.  Strength  cr„ 

-25.80 

-20.89 

ai 

(m) 

h 

(mm) 

HQQIIIIII 

71 

2.2 

2.2 

2.8 

Roughness  Parameters 

RMS  roughness  Ck 

(cm) 

0.87 

0.96 

(0.96) 

(0.96) 

Spectral  strength  ^ 

(cm^) 

0.00081 

0.00175 

(0.00175) 

(0.00175) 

Spectral  exponent  7 

3.02 

3.21 

(3.21) 

(3.21) 

Gradient  Parameters 

Density  po 

(g/cm^) 

1.52 

1.37 

1.18 

1.1 

Sound  speed  cq 

(m/s) 

1550 

1510 

1495 

1515 

a 

(m“*) 

300 

800 

800 

100 
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4 

Model-Data  comparison 


Acoustic  data  were  acquired  using  a  single  element  transmitter  and  a  multi-element 
receiver  mounted  on  a  pan/tilt  mechanism  installed  at  the  top  of  a  telescopic  tower. 
The  source  height  above  the  seabed  was  approximately  4  m.  An  initial  survey 
using  a  300  kHz  multibeam  bathymetric  echosounder  ensured  that  the  large  scale 
slope  was  negligible.  Sidescan  sonar  data  and  diver’s  observations  confirmed  the 
spatial  stationarity  of  the  superficial  seabed  properties  over  the  acoustic  footprints. 
Reception  on  the  64  element  receive  array  enabled  an  angular  beamforming.  The 
backscattering  strength  as  a  function  of  grazing  angle  was  computed  by  averaging 
independent  beams  over  a  receive  sector  width  of  30°.  Data  displayed  below  are  for 
a  transmission  of  50//S  CW  pulses  with  a  center  frequency  of  140  kHz). 


4,1  Effect  of  the  gradient  on  the  reflection  coefticlent 

Figure  24  illustrates  the  effect  of  vertical  gradients  on  the  reflection  coefficient  at  each 
of  the  four  experiment  sites.  Near  normal  incidence,  the  effective  reflection  coefficient 
is  lowered  and  approaches  the  reflection  coefficient  obtained  for  a  sediment  having  an 
impedance  of  pqcq  whereas  at  low  grazing  angle  and  particularly  around  the  critical 
angle,  the  effective  reflection  coefficient  gets  closer  to  the  reflection  coefficient  for  an 
impedance  of  pgCs. 


4.2  Tellaro 

At  Tellaro,  as  shown  in  Fig.  25,  for  the  case  of  no  gradient,  modeled  interface 
backscattering  strengths  and  volume  backscattering  strengths  are  respectively  too 
high  and  too  low  compared  to  the  measured  backscattering  strengths.  On  the  con¬ 
trary,  in  the  presence  of  property  gradients,  as  the  reflection  coefficient  is  lowered 
by  about  3dJ5,  the  modeled  interface  backscattering  strength  matches  the  data  well. 
However,  scattering  within  the  volume  is  clearly  not  a  competing  mechanism  at  this 
site  and  at  this  frequency.  This  can  be  explained  by  the  low  amount  of  scatters 
within  the  volume  and  by  a  relative  higher  impedance  contrast  compared  to  the 
other  sediments  studied  in  this  paper. 
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4.3  Venere  Azzura 

At  Venere  Azzura,  interface  and  volume  scattering  have  comparable  levels  (see 
Fig.  26).  In  the  absence  of  gradients,  interface  scattering  dominates  at  all  graz¬ 
ing  angles  and  is  slightly  higher  than  the  measured  level.  The  sum  of  interface  and 
volume  scattering  overestimates  backscattering.  However,  if  vertical  gradients  are 
taken  into  account,  interface  scattering  is  reduced  by  about  SdB.  The  sum  of  in¬ 
terface  and  volume  scattering  matches  almost  exactly  the  measured  backscattering 
strengths.  In  the  low  grazing  angle  portion,  the  measured  backscattering  strength 
decreases  rapidly  as  the  grazing  angle  decreases.  This  effect  is  not  observable  in  the 
model  even  though  self-shadowing  is  taken  into  account.  One  may  actually  question 
the  physical  processes  responsible  for  the  decay  of  backscattering  strength  at  low 
grazing  angle  as  well  as  accuracy  of  the  simple  model  considered  here.  Clearly,  self¬ 
shadowing  is  partly  responsible  for  this  decay  even  though  roughness  is  often  low 
for  fine  sediments,  but  it  docs  not  seem  to  be  the  only  cause.  As  a  suggested  com¬ 
peting  mechanism,  the  spatial  fluctuation  of  sound  speed  at  the  sediment  interface 
could  also  be  responsible  for  this  rapid  decay,  especially  for  highly  inhomogeneous 
sediments.  The  presence  of  a  distribution  of  sound  speeds  affects  local  refraction,  re¬ 
flection  and  transmission  with  a  subsequent  impact  on  the  backscattering  strength. 
The  effect  of  local  interface  sound  speed  fluctuation  is  not  treated  here  but  should 
be  the  object  of  further  investigations. 


4.4  Porto  Venere 

Backscattering  at  Porto  Venere  is  high  and  higher  than  at  TE  and  VA,  despite  a 
similar  low  roughness  and  a  lower  water-sediment  impedance  contrast.  This  level  can 
only  be  explained  by  the  higher  velocity  and  density  fluctuation  within  the  volume. 
As  can  be  seen  in  Fig.  27,  volume  scattering  dominates  at  moderate  and  low  grazing 
angles  both  in  the  presence  and  in  the  absence  of  gradients.  Given  the  low  water- 
sediment  density  contrast  the  effect  of  the  vertical  gradient  is  even  more  significant 
in  reducing  the  relative  interface  scattering.  In  the  prevscncc  of  property  gradients, 
volume  scattering  is  unambiguously  the  only  significant  scattering  mechanism. 


4.5  Punta  della  Mariella 

The  physical  mechanisms  responsible  for  scattering  at  Punta  della  Mariella  are  the 
same  as  at  PV:  sediments  are  strongly  inhomogeneous  and  have  a  low  impedance.  In 
the  absence  of  gradients,  the  sum  of  interface  and  volume  backscattering  strengths 
over  predicts  the  measured  levels.  However,  with  a  gradient,  interface  scattering  is 
reduced  and  volume  scattering  alone  is  very  close  to  the  data.  The  slight  overprcdic- 
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tion  of  volume  scattering  alone  may  be  due  to  an  inaccurate  estimation  of  some  of  the 
volume  geoacoustic  properties  (gradient,  inhomogeneity  levels,  degree  of  correlation 
between  density  and  velocity  fluctuation).  Internal  seabed  properties  are  retrieved 
from  a  small  volume  (10x10x36  cm^  at  PdM)  which  is  assumed  (possibly  in  error) 
representative  of  the  properties  of  the  entire  acoustic  footprint.  However,  general 
agreement  between  the  model  including  gradients  and  the  acoustic  backscatter  data 
is  excellent  and  consistent  for  the  4  sites  considered. 
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Telloro 


Venere  Azzura 


Punto  Dello  Moriella 


Figure  24  Effect  of  the  vertical  gradient  on  the  reflection  coefficient.  The  effective 
reflection  coefficient  {c{z)^  p{z))  is  com,puted  using  OASR. 
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Venere  Azzura  without  gradient 


Figure  26  Model- Data  comparison  for  the  Venere  Azzura  case. 
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Porto  Venere  without  gradient 


Figure  27  Model-Data  comparison  for  the  Porto  Venere  case. 
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Punta  Della  Moriella  without  gradient 


Punta  Della  Moriella  with  gradient 


Figure  28  Model-Data  comparison  for  the  Punta  Della  Mariella  case. 
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5 

Summary 


Four  sites  in  the  Gulf  of  La  Spezia  (Italy)  have  been  carefully  geoacoustically  char¬ 
acterized  using  in  particular  X-Ray  computed  tomography  on  cores  and  digital 
photogrammetry  on  stereo  photograph  pairs.  This  thorough  characterization  has 
allowed  a  simple  backscattering  model  to  be  constrained.  A  comparison  between 
backscattering  data  versus  angle  acquired  at  four  sites  comprised  of  unconsolidated 
sediments  and  the  developed  model  has  been  presented.  For  two  of  the  sites  (PV 
and  PdM),  volume  backscattering  from  velocity  and  density  fluctuations  explains 
the  high  backscattering  levels.  It  also  suggests  that  consideration  of  gradients  in 
density  and  sound  speed  are  necessary  to  explain  the  measured  levels.  Regarding 
the  correct  way  of  treating  scattering  from  an  inhomogeneous  medium  bounded  by 
a  rough  interface  and  affected  by  gradients,  many  questions  still  remain  to  be  an¬ 
swered.  It  is  also  evident  that  at  high  frequency,  in  situ  geoacoustic  measurements 
of  high  spatial  resolution  and  accuracy  become  more  delicate  and  important. 

We  can  conclude  with  the  following  comments: 

•  The  density  structure  of  highly  bioturbated  sediments  does  not  seem  to  behave 
like  a  modified  power  law  in  the  low  spatial  frequency  portion  of  the  power 
spectrum.  However,  at  high  spatial  frequencies  (typically  above  logio{k)  >  2.8 
which  corresponds  to  acoustic  frequencies  higher  than  75kHz  in  the  backscat¬ 
tering  case,  the  power  spectrum  decay  can  be  acceptably  described  by  a  power- 
law.  At  high  spatial  frequencies,  the  power  exponent  of  the  density  structure 
71  seems  to  increase  as  the  level  of  fluctuation  increases:  (TE:  71  =  2.2,  VA: 
71  =  2.2,  PV:  71  =  2.8,  PdM:  71  =  3.0).  However,  the  size  and  number  of  the 
core  samples  may  not  be  sufficient  to  draw  conclusions  on  the  general  behav¬ 
ior  of  the  density  power  spectrum  of  bioturbated  sediments,  especially  at  the 
lower  frequency  end.  More  and  larger  core  samples  are  needed  to  draw  general 
conclusions  at  lower  spatial  frequencies. 

•  At  the  most  bioturbated/het erogenous  sites  (PV  and  PdM)  and  in  the  mod¬ 
erate  to  low  grazing  angle  region,  the  dominance  of  scattering  is  essentially 
caused  by  a  high  level  of  internal  geoacoustic  property  fluctuation. 

•  In  this  study,  the  sediment  density  structure  of  highly  bioturbated  fine  sedi¬ 
ments  is  mainly  isotropic  in  the  upper  10  cm. 
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•  Vertical  gradients  have  a  major  effect  on  backscattering,  especially  as  the 
acoustic  frequency  increases.  When  the  acoustic  wavelength  is  of  the  order  or 
smaller  than  the  vertical  size  of  the  gradient,  the  relative  contribution  for  the 
rough  interface  to  the  total  scattering  process  is  significantly  reduced.  In  the 
absence  of  gradients,  model-data  comparisons  arc  not  satisfactory.  However, 
when  gradients  arc  taken  into  account,  model- data  comparisons  arc  greatly 
improved. 

•  The  simple  models  presented  here,  despite  the  general  agreement  with  the 
measured  backscattering  strength  may  not  be  accurate  in  the  very  low  grazing 
angle  region.  They  do  not  account  for  mechanisms  such  as  anisotropy  of  the  in¬ 
terface,  interface  and  volume  multiple  scattering,  the  coupling  between  a  rough 
boundary  and  an  inhomogenous  medium,  combined  density /velocity /attenuation 
gradient  effect,  bioturbation,  discrete  scatters,  gradients  in  the  heterogeneity 
statistics,  etc.  These  aspects  warrant  particular  attention  in  future  research 
works. 

•  It  should  also  be  stressed  that  even  though  current  technology  (e.g.  digital 
photogrammetry  and  X-Ray  tomography)  enables  a  better  statistical  descrip¬ 
tion  of  the  geoacoustic  environment,  this  technology  may  still  remain  insuffi¬ 
cient  above  100  kHz  given  the  complexity  of  the  considered  seabeds  and  the 
scattering  mechanisms  taking  place  in  them. 
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